Representation of entanglement by negative quasi-probabilities 
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Any bipartite quantum state has quasi-probability representations in terms of separable states. 
For entangled states these quasi-probabilities necessarily exhibit negativities. Based on the general 
structure of composite quantum states, one may reconstruct such quasi-propabilities from exper- 
imental data. Because of ambiguity, the quasi-probabilities obtained by the bare reconstruction 
are insufficient to identify entanglement. An optimization procedure is introduced to derive quasi- 
probabilities with a minimal amount of negativity. Negativities of optimized quasi-probabilities 
unambiguously prove entanglement, their positivity proves separability. 
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I. INTRODUCTION 



Among the most striking discrepancies between clas- 
sical and quantum physics are the nonclassical correla- 
tions between the subsystems of a compound system [l[ . 
In general, a common description of all subsystems is 
needed for describing the properties and the evolution 
of the compound system [2| ■ These nonclassical correla- 
tions are usually called entanglement. Today, entangle- 
ment is considered to be the key resource of the vast fields 
of Quantum Information Processing, Quantum Compu- 
tation, and Quantum Technology, see e.g. [1, Q|. The 
most prominent examples of entangled states are the 
Bell states, which maximally violate Bell-type inequal- 
ities which are valid for classical correlations [f|. 

On the other hand, in the field of Quantum Optics non- 
classical effects have been studied over many years, which 
are related to negativities in the Glauber-Sudarshan dis- 
tribution H, 0], in the following denoted as Pqs- It de- 
scribes any quantum state of the harmonic oscillator as 
a pseudo- mixture of coherent states. A state is consid- 
ered to be classical, if the Pes function has the properties 
of a classical propability distribution. If this is not the 
case, then the state is called a nonclassical one @, Q. 
Recently, negativities of Pes could be experimentally 
demonstrated 

Entanglement and the nonclassical effects reflected by 
the negativities in Pqs appear to be related in some form. 
When mixing two radiation fields by a beam splitter, 
nonclassicality in one of the input channel is necessary 
for the gen eration of entanglement in the output chan- 
nels [ITI Il2| . On the other hand, even an entangled com- 
posite state has negativities in the Pes function, which is 
a necessary but not sufficient condition for entanglement. 

A detailed understanding of the structure of entangled 
quantum states is a topic of fundamental interest. The 
idea of the decomposition of entangled states into a best 
approximation by a separable state and a remaining cn- 
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tangled state delivered further insights into the nature of 
entanglement, see did. It has been shown in [H Eft 
that any quantum state can be given in terms of local 
pseudo-mixtures, here and in the following denoted by 
the quasi-probability distribution PEnt, that visualizes 
entanglement. An interesting consequence is in this con- 
text, that local measurements are sufficient to describe 
any quantum system completely (TtI [Tc| . including en- 
tanglement. 

In the present contribution we reconsider the represen- 
tation of entangled states by quasi-probabilities, PEnt, 
in terms of separable states. This representation is 
compared with the Glauber-Sudarshan representation of 
nonclassical states by classical ones. We give an ex- 
plicit construction scheme for Pent for arbitrary quantum 
states, which opens the possibility of reconstructing these 
quasi-probabilities from experimental data. However, the 
quasi-distributions resulting from the bare reconstruction 
are not unique. The identification of entanglement would 
require to inspect an infinite set of quasi-probabilities. 
To overcome this problem is the main aim of our pa- 
per. Starting from the problem that PEnt is not unique, 
we derive an optimization procedure to obtain PEnt, opt, 
representing a quasi-probability of a given state whose 
negativities are properly minimized. The most impor- 
tant consequence is that a quantum state is entangled, 
if and only if PEnt,opt has negativities. The method can 
be applied to arbitrary quantum states. It identifies if 
a state is separable or entangled by a single optimized 
quasi-distribution. 

The paper is structured as follows. In Sec. |TT] we re- 
consider both the description of nonclassicality by neg- 
ativities in the Glauber-Sudarshan representation and 
the corresponding description of entanglement by neg- 
ativities in the related quasi-probabilities. The problem 
of the reconstruction of quasi-probabilities characteriz- 
ing entanglement from measured data is considered in 
Sec. IIIII In Sec. IIVI we derive an optimization method 
for the quasi-probabilities, which is based on the separa- 
bility eigenvalue problem [19| . We apply this approach to 
a separable and an entangled quantum state. A summary 
and some conclusions are given in Sec. fVl 
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II. QUASI-PROB ABILITY REPRESENTATIONS 

Let us start with the well known Glauber-Sudarshan 
representation of a single-mode quantum state as an 
example for quasi-probabilitics in statistical quantum 
physics. A coherent state can be introduced as the eigen- 
state of the annihilation operator a of a harmonic oscil- 
lator, a\a) = a\a) (a is a complex number). It has been 
shown in 0, 0] , that any state can be given in the form 
of the Glauber-Sudarshan representation 



(1) 



P = dP GS (a)\a)(a\, 



with Pqs being a signed measure, or more precisely a 
quasi-propability distribution. By definition a state is 
nonclassical, if Pes fails to be represented as a classical 
propability distribution P c i ass , see Thus, a nonclas- 

sical state p has negativities in its Pqs distribution. The 
only classical pure state (with a non-negative distribution 
Pqs) is the coherent state [2(|. Thus the implication of 
Eq. (p} is, that any quantum state p can be given solely 
in terms of classical states |a). 

As mentioned above it has been shown in [lfj [l(| , that 
any quantum state for finite Hilbert spaces can be given 



p = (l + p)a-pa', 

with <T and a' being separable states and p > 0. 
condition for positivity, (ip\p\ip) > 0, leads to 



V|V) : (V>|ff|V> > 



1 



/' 



■W\i>)- 



(2) 
The 



(3) 



A state p is separable, if it can be written as p = a. This 
means p = 0. There are two items describing the failure 
of an entangled state g to be a convex combination of 
separable states. For the latter it is not allowed to use 
coefficient with both (1 + p ) > 1 and (—pi) < 0. 

As we have shown in |2l|, only finite systems need to 
be considered to characterize entanglement completely. 
Hence, it is not a restriction, if one assumes a finite, but 
arbitrary dimensional Hilbert space H.a'&'Hb- Thus, any 
quantum state can be given as 



p = ^Pk\ak,b k )(ak,bk\, 

k 



with the quasi-probabilities 

PEnt = (Pk)k- 



(4) 



(5) 



The negative quasi-propabilities can be considered to be 
the key signature of entanglement. 



III. RECONSTRUCTION OF 
QUASI-PROBABILITIES 

This leads us to the question, how to obtain the quasi- 
probability distribution PEnt ~~ in particular for entan- 
gled states. We will start to consider two quantum sys- 
tems represented by Ha and H.b, given by the bases 
{|fc)}fe=i,2... for each system. A quantum state can be 
given, for example, as \a, b) = \a) ® \b). The states |o) 
and | b) describe the system A and P, respectively. The 
linear structure of quantum theory allows us to superim- 
pose quantum states. Let us consider the superposition 
state 



1 



V2 



(M + IM)). 



(0) 



Obviously this state is not factorizable. 

To show how to describe entanglement with quasi- 
propabilities, we start with the Bell state g = |$ + ) (<f> + 1. 
The statistical operator g corresponding to this state is 
given as 



g =|$+)($+ 
1 



(7) 



^ (|0,0)(0,0| + |1,1)(1, 1| + |0,0)(1, 1| + |1,1)(0,0|) . 

To rewrite this state in terms of separable states, we need 
to define for each system the states \s n ), 

| S „) = -L(|0)+H1)) (n = 0,l,2,3), (8) 

and for these states we define factorizable statistical op- 
erators <3„ = |s n )(s„| <S> \s n )(s n \. Now we can obtain a 
superposition, operator €, of these separable states 6 n 
with the real coefficients ±1, 



<£ = 6 - 6i + 6 2 - e 3 
= |0,0)<1,1| + |1,1)<0,0|. 



(9) 



Therefore, we can rewrite the density operator for the 
Bell state |$ + ) in terms of separable states, 



e=|$+)(a>+| = -(|o,o)(o,o| + |i,i)(i,i| + e) 

=-|0,0)(0,0| + -|1,1)(1, 1| + i|s ,So)(so,s | - i|si,s 1 )(si,s 1 | + | |s 2 , s 2 ) <s 2 , s 2 1 - -|« 3j *3><*3.*3|. (10) 
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For the Bell state |5> + ), from Eq. (JTUJ) we see that the 
quasi-propability distribution PEnt has negativities, 



(p 



k)k 



1 1 1 

2' 2' 2' 



1 1 

2' 2' 



(11) 



This is a simple example of a quasi-probability represen- 
tation of an entangled state by separable ones. 

Now we will show that the simple example of a Bell 
state is sufficient to reconstruct any mixed quantum 
state. We may rewrite the matrix pkimn, representing 
a general quantum state, in terms of separable quantum 
states \a k , bk){a.k, b k \ and a quasi-propability distribution 
■Ffent — (Pk)k- Let us consider a density operator of a bi- 
partite system, 



P = Pkimn\kJ)(m,n\ 

klmn 



(12) 



that has been obtained from experimental data. By solv- 
ing the eigenvalue problem of the matrix p, we obtain the 
spectral decomposition 

P = ^2p^\4>i}(&\, (13) 

i 

with p r t )i > and J^p^ = 1) which is a common proce- 
dure. 

It is sufficient to show how to write any pure state 
\ip)(ip\ solely in terms of separable states. For any pure 
state one can obtain the Schmidt decomposition [3J] 



£ 

k=l 



Afe \e k ,f k ) 



(14) 



with Xk > 0, the Schmidt rank r(ip), and {|efe)}fe and 
{\.fk')}k' being orthonormal in Ha and He, respectively. 
Under the assumption \ip) ^ \a) ® |&), we can write 



= S A *Ai|e fcj / fc )(e i ,/,| 

= 5>S|e fcl / fc >(e fc ,/ fc | ( 15 ) 

k 

+ ^A fc A ; {\e ki f k )(ei,fi\ + \e u fi)(e k) f k \) , 



The first sum, J^ fe A||efc, fk)(sk, fk\, is a separable state. 
The second sum depends on the terms (for k > I) 



&k,l =\ek,fk){ei,fl\ + h, fi)(e k Jk\ 



|4^><4' 



(k,i). 



|S (M) )(S (M) 



(16) 



|4 M) ><* 



(*.0 1 

2 



with |sV ) = 7f d e fe' /fc) + *"l e '> /'))■ lt can be written 
analogously to € as a superposition of separable states 
\sit' 1 ^) with positive and negative coefficients. According 
to Eq. , we can give an expansion of any pure state in 
terms of separable states, IV'KV'I = (1 + /v)^ — MV" 1 ^' 
with jtiy, = 2^] fc >; A fc A ! an d 



a i> -TT I 2^ A k\ e k, Jk){ek, Jk\ + 2^ A fc A ' (J S ' S A s > S I + l S 2 ' S 2 ) ( S 2 > S 2 IJh ( 17 j 

^ \ k k>l / 

+ i4 M) ,4 M) )(4 W) ,4 W) l) • (is) 

I 



By using Eq. (|13| . any mixed quantum state p can be 
written as a convex combination of pure states. There- 
fore, we conclude that any state p can be written as 

i i 

= J2lk\ak,b k )(ak,b k \, (20) 
k 

solely in terms of separable states. 

We have rewritten the matrix p k i mn , that may be 
reconstructed from experimental data, in terms of the 



separable quantum states \a k , bk}(ak, b k \ and the quasi- 
propability distribution PEnt = {q k )k- In principle, one 
may obtain the latter by the procedure presented above. 
By experimental state reconstruction one gets pkimn , the 
spectral decomposition yields p ( j >i and \4>i), via Schmidt 
decomposition one obtains Xj and \ej,fj), which finally 
leads to the quasi-probabilities qk and \a k ,b k )- 

An arbitrary quantum state can be expressed in an 
integral form, 

P= f dP Ent (a,b)\a,b)(a,b\, (21) 
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by the non-orthogonal states \a, b)(a, b\ and the quasi- 
propabilities distribution PEnt characterizing entangle- 
ment. Negative values of PEnt are necessary signatures of 
entanglement. The form of PEnt - as a signed measure - 
can be written by the reconstruction scheme given above 
as 

Pnnt = YllkS\ ak ,b k ), (22) 

k 

with the Dirac measure $\ ak ,b k )- 

The integral form of the state p in terms of PEnt , see 
Eq. pTj) . is analogous to the Glaubcr-Sudarshan rep- 
resentation with coherent states, see Eq. (fT]). In both 
cases the failure of a classical interpretation of PEnt and 
Pgs describes the quantumness of a state. Let us as- 
sume we have a system of two harmonic oscillators and 
PEnt (a, b) — for all factorizeable states \a,b) which are 
not coherent ones, \a,b) ^ \a,(3). In this case we obtain 
PEnt = Pgs, with negativities for entanglement. Further- 
more it is obvious that negativities of Pqs are not neces- 
sarily due to entanglement. For example, let us consider 
the Fock state \n A ) <£> \ns) (nA,riB > 1), 

p=\n A ){nA\®\n B ){nB\- (23) 

This is a nonclassical state, but obviously a separable 
one, P En t = Pciass = &\n A )&\n B )- We conclude that entan- 
glement is a subclass of all nonclassical phenomena, see 

Kg. in 




FIG. 1: Relation between entanglement and nonclassicality of 
quantum states. A negative Pgs quasi-distribution is neces- 
sary for entanglement (light gray area). A classical Pgs dis- 
tribution is sufficient for separability (dark gray area). There 
are quantum states which are separable but not classical ones, 
see Eq. J23J. 



The decomposition of a quantum state as a convex su- 
perposition of and \4>){4>\ is not unique. Thus, the 
distribution pEnt is also not unique, but at least one re- 
construction scheme has been presented in this section. 



Let us denote by Af the set of all signed measures /, 
generating the operator (for all \if>): 0\ip) =0), 

6= / df(a,b)\a,b)(a,b\. (24) 

For all / and /' in M and all real numbers v and v' , 
the signed measure vj + v' f generates again 0. We can 
write p = p + 0, or we obtain a new measure PEnt + / 
for any element / of Af, that generates the same state. 
Under all possible PEnt, which generate a given state p, 
the quasi-propability PEnt, opt with some kind of minimal 
amount of negativities needs to be found, 

j d(\P Ent (a, b) + g{a, b)\ 2 ) min . (25) 

This requires an optimization as presented below in 
Sec. IIV1 Obviously, any possible quasi-probability Ps n t 
for an entangled state g has negativities and a state g is 
entangled, if and only if PEnt fails to be a classical propa- 
bility distribution for any decomposition, PEnt 7^ Pciass ■ 
Therefore entanglement is due to negativities in the dis- 
tribution PEnt,opt- 

The other way around, as it was known for two 
decades (22|, any separable state can be given by a clas- 
sical distribution, Psnt = Pciass- Thus a single non- 
negative quasi-distribution PEnt is sufficient to demon- 
strate the entanglement of a given quantum state. The 
advantage of an optimized quasi-probability PEnt, opt is 
that it clearly uncovers entanglement and separability 
by its negativities and its non-negativity, respectively. 



IV. OPTIMIZATION 
A. Derivation of Optimization 

The Hilbert space Hab = H A ® H B = C dA <g> C dB 
is finite dimensional, dimH ab — d A d,B- As we have 
mentioned above, the assumption that our Hilbert space 
is finite dimensional is not a restriction for the property 
of entanglement [2l|. The Banach space Herm(7i J 4s) of 
Hermitian Operators for the given Hilbert space is a finite 
dimensional real vector space, dim}ierm(li. ab) = Aa^b- 
First of all we want to introduce the separability norm 

II ' 1 1 SCp 5 

||i|| SC p = Bup{Ka,6|A|o,6)| : (a,b\a,b) = 1}. (26) 

Now we prove that the definition given in Eq. (|26p is a 
norm for Herm^^s)- Obviously ||A-A|| scp = |A|||A|| sep 
and ||i + P|| scp < ||i|| scp + ||P|| scp hold true, 

(a, b\\A + nB\a, b) = X(a, b\A\a, b) + n(a, b\B\a, b) 
|(o,6|xA + /6B|o I 6)| < |A||(a,6|i|a,6)| + |k||(o,6|B|o,6)| 

||Ai + KP||scp < |A|||i|| S cp + |K|||P|| se p. 
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Further on ||j4|| sep = is equivalent to A = 0, 
||i||sc P = & y\a,b) : (a,b\A\a,b) =0 
VP E „t : J dP Ent (a,b)(a,b\A\a,b) =0 

& V|V> : (il>\A\i>) = & A = 0. ■ 

Let .M be the set of separable quantum states. This 
set is bounded and closed. The set A^o = {\a,b)(a,b\ : 
(a,b\a,b) = 1} is the set of pure separable states. The 
set of all separable quantum states can be defined as the 
set of all possible convex combinations of pure separable 
states, M = convex (.Mo)- 

A quantum state p is separable, if and only if it exists 
a separable quantum state a such that the distance of p 
and this state is zero, 

3<r = ^2pk\xk,Vk)(xk,yk\ ■ 

k 

Pk>0, \\p-a\\ sep = 0. (27) 

The norm || • || sop is the least upper bound (sup) of the 
projections (a,b\ ■ \a,b). Thus, Eq. (f27|) is equivalent to 
the optimization 

g(a,b) = (a, b\p- &\a,b) -> g opt (28) 

under the condition h(a,b) = (a,b\a,b) = 1 and the as- 
sumption that all optimal (especially the largest and the 
smallest) values are g op t = 0. We discussed an optimiza- 
tion procedure as it is given in Eq. ([28]) in For this 
purpose we derived a set of equations - the separability 
eigenvalue equations - for a linear operator i, 

L b \a) = g\a) 

L a \b)=g\b), (29) 

with (a,b\a,b) = 1, L b = ti B (L [i A (g> \b)(b\]) and L a = 

ti A (L [|o)(o|®i B ]). 

Obviously the norm is zero, \\p — er|| S cp = 0, if and 
only if all solutions (<?„, \a n , b n }) of the separability eigen- 
value problem for p are solutions for a and the other way 
around, 

Pb n \a n ) = g n \a n ) 

Pajb n ) = g n \b n ) (30) 
<=> 

<Tb n \<hi) = gn\a n ) 

&ajb n ) = g n \b n ). (31) 

Thus we have to solve the separability eigenvalue problem 
of the Hcrmitian operator p = a. We obtain the following 
equations: 

g n = (a n ,b n \a\a n ,b n ) = Y^Pfe \(a„,b n \x k ,y k )\ 2 , (32) 

k G„, k 



with g = (g n ) n , P = (Pk)k and G = {G, hk ) n ,k this is 
equivalent to 

g = Gp. (33) 

At least one solution pb = {pk)k — PEnt must exist for 
any p, since any state can be written in terms of separable 
states as presented in Sec. Mil The state p is separable, 
if and only if 

3p solution of Eq. ([33]) : (p)„ = p n > Vn. (34) 

The kernel T = ker(G) gives us the possibility to obtain 
such a positive vector. It is related to the measure / 
given in Eq (|24|). / G Af. Let 7 G T be an arbitrary 
element of the kernel. Thus po + 7 is another solution 
of Eq. (|33[) . If we rewrite this in therms of measures we 
obtain / = 7. The quantum state a is separable, if and 
only if (po + t*)„ > for all n for one choice for 7. 

The question arises how to obtain the matrix elements 
G n ,k = \(a n ,b n \x k ,yk)\ 2 or what are the states \x k ,y k ). 
If we would know the linear map G we could identify 
its kernel T and all solutions p. This could solve the 
separability problem for a given quantum state p. 

For pure separable quantum states a = \x, y) (x, y\ the 
only non-trivial solution of the separability eigenvalue 
problem is (l,\x,y)). Thus \x k ,y k ) = \x,y) for g k = 1 
fullfills the desired conditions. Ignoring this trivial case, 
let us assume a is a mixed and separable quantum state 
Afo). 

As we have mentioned above, the set of separable quan- 
tum states M. is a convex set. Thus it is a subset of a 
hyperplane of Herm(7i ab)- Let us shift this convex set 
M! = M — a. Thus the state under study is &' = 0. The 
subspace of Herm(H J 4s) given by V = span(A4') and a 
norm restricted to this space, || • ||y', is a normed vector 
space. The boundary of M! with respect to the vector 
space V" is dM! = A4o — &. Thus M! has a finite volume 
not equal to zero, M'\dM' ^ 0. 

Now it is usefull to find the extremal elements of the 
convex set dM' for the shifted state under study, a' = 0, 

ll^purellv -»■ optimum (cr purc G dM'). (35) 

Let us denote this optimal elements by (7 pure i ■ Since the 
M' has a non-zero and finite volume, there are elements 
i k such that all cr purc 4 are linearly independent (for ex- 
ample in Fig. [2] the vectors from the state to the gray 
dots 1 and 3), and other elements ii such that p' puro ^ il is 
linearly independent (in Fig. [5] the vectors from the state 
to the gray dots 2 and 4). Further on these elements for 
i k and ii give boundaries for M' in all needed directions. 
This means that the tangential hyperplanes of dM' at 
the points given by i k and ii describe a closed set around 
M '. Since all &' puietik and all <j' purcA[ are directed in this 
way, the point a' = can be given as a convex combina- 
tion of all a' pmc tk and o- pm . GA , 

G convex({<T purc lk } k U {a' pweM }i). (36) 
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With other word we can say, that all extremal points 
of a' at the boundary dM.' generate a convex set 
convex({<Tp Uro which includes this clement a', see 
Fig.H 




1 



• separable quantum state 

• separability eigenvector 



FIG. 2: The light gray area is the set of separable quantum 
states M . The radii of the dashed circles are the separability 
eigenvalues. The dark gray area is the convex set given by 
the separability eigenvectors. The separable state is element 
of this set. 



Since Herm(H.AB) is finite dimensional, it follows that 
the norm || • || scp is equivalent to any other norm, for 
example given by || ■ ||y. Thus the optimization can 
be performed for the separability norm, which leads to 
the separability eigenvalues. Now we can undo the shift- 
ing. Therefore the solutions of the separability eigen- 
value problem for a yields our desired vectors |afc, bfc) = 
\%kiDk)- Thus our map G is given by 



Gk,i = \{a k ,b k \ai,bi)\ 2 

= tr (\a k ,bk)(ak,b k \\at,bi)(ai,bi\). 



(37) 



The map G is symmetric and positive semi-definite, since 
for all real vectors if holds 



t t Gtt 



k,l 



T^kGkA^l 



=tr 



=tr 



^ TTk\ak,b k )(a k ,bk\ 

k 

y"Vfc|afc,frfc)(a fc ,bfc| 



y^^TTi\ai,bi)(ai,bi\ 



> 0. 



(38) 



Therefore all eigenvalues of G are positive or zero, and 
a spectral decomposition of G with real eigenvectors is 
possible. 



In conclusion we showed that the matrix G is given 
by the separability eigenvectors of the state p = a. The 
matrix G itself is a positive semidefinite symmetric oper- 
ator and therefore all eigenvectors and eigenvalues can be 
calculated with algorithms known from linear algebra. If 
our state under study is not a convex combination of pure 
factorizeable states (inseparable), the procedure fails to 
deliver a solution with (p) n > for all n. This is due to 
the fact, that the state is not element of the convex set 
of separable states A4. 



B. The Optimization Procedure 

Now we want to describe the optimization scheme as 
derived above. We find PEnt, opt for a given quantum state 
p in the following way. 

1. We solve the separability eigenvalue problem of p, 

p b \a) = g\a) 
Pa\b)=g\b). 

The solutions are given by (g n , \a n ,b n ))- 

2. The vector g = (g n ) n and the linear map G = 
(|(a/c, bk\ai, bi)\ 2 )k,i can be calculated. Now the (ordi- 
nary) eigenvalue problem of G can be solved, 

Gpx = ApA- 

A solution of Eq. (|33|) . g = Gp, can be found. The result- 
ing quasi-distribution for entanglement is PEnt — (Pn)n = 
p, and 



}(a n ,b n \ 



3. The orthogonal eigenvectors p Ufc for the eigenvalue 
A = are the basis for the kernel T. Thus the optimized 
-PEnt, opt — Popt and it can be given as 



Popt = P 
with real coefficients Cfc, 

c k = - 



k 



CkP0 k , 



Po k P 

P^ k PO k 



The optimized quasi-distribution is PEnt, opt — Popt- 

Let us prove the given form for Ck- According to 
Eq. (f25|) . the minimal amount of negativities needs to 
be found, 

^2 KPopOrJ 2 -> min. 

n 

This is equivalent to the optimization for (ck)k, 

C((c k )k) = ^2(Popt)l =PoptPopt 

n 

CkPo k \P 



C kP0k 
k 
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The optimum can be calculated for each q as 
dC 



da 



C. Examples for the Optimization 

To illustrate the strength of our method, we consider 
two examples: a mixed separable state and a Bell state. 
For obtaining analytical solutions, we choose simple ex- 
amples. In General, the given optimization procedure 
can be implemented with numerical methods. However, 
the examples provide some basic properties of the opti- 
mization. 

Example A. The first state a is separable, it is a con- 
vex combination of five pure separable states, 

<T=i|0,0)(0,0| + i|l,0)<I,0| + i|0,I)(0,I| 



glMXMl + ^|so,So>(so,s | 



(39) 



with |so) as defined in Eq. {5}. To solve the separability 
eigenvalue problem of a, we use the Proposition 1 given 
in 0, 

a\a,b) = ^\a,b) + ^(a + ai)(&o + h)\s Q , s ), (40) 

with | a) = a |0) + aijl) and \b) = b \0) + h\l). The 
vector \a,b) is a separability eigenvector, if a\a,b) has a 
Schmidt decomposition with \a, b). We can conclude that 
the solutions (<?, \a, b)) are 



77; l s 2 , S ^ 



77, l s 0,S2j 



Now we obtain <?, G and the solution of Eq. (|33l) , 



1 



,g=-(5,l,l,l) J 



G = 



/ 1 

10 

10 

\ 1 



(41) 



i>= 1(5,1,1 A) T = P< 



cnt,opt- 

The optimization map G has no eigenvalue equal to zero. 
Thus an optimization becomes superfluous. We conclude 
that the optimal decomposition of a consists of only four 
separable pure states, 

5 1 

a =-|s ,s )(so,So| + g -|s 2 ,So)(s 2 ,So| 

+ Tj|S0,S2)(S0,S2| + 7jl s 2,52)(S2,S 2 |. (42) 



This example also shows, that the initial pure separable 
states and their coefficients need not necessarily be the 
same as the optimization procedure delivers. 

Example B. After this simple example let us consider 
a more sophisticated one. Therefore let us consider the 
previously studied Bell state g = |$ + )(i> + |. Its impor- 
tance is due to the reconstruction as presented in Sec. lIII| 
which depends only on the reconstruction of this Bell 
state. 

For the solution of the separability eigenvalue problem 
of q, let us perform the same procedure as presented in 
the first example, 

g\a, b) = ^(a Q b + ai&i)(|0, 0) + |1, 1)). (43) 
We obtain the non-trivial solutions, 



|,Ia>® |a*> 



and the trivial ones satisfy 

a b + ai&i = 0. 



(45) 



If the state is separable, PEnt.opt (a, b) = F c i ass (a, b) > 
0, then we can conclude that all needed separabil- 
ity eigenvectors are orthogonal to all \a,b) satisfying 
Eq. ([45]) . For the example \a,b) = |0, 1) we obtain 



(0,l|g|0,l)= / dP class (a)(0,l\(\a,a*){a,a*\)\0,l) 

rfPcia SS (a)|(0|a)| 2 |(a*|l)| 2 = 0. (46) 



This implies P c i ass (a) = or (0|a) = or (l|a*> = 0. 
This means the only needed separability eigenvectors are 
|0,0) and |1,1), but 

^ ( 7=i|0,0)(0,0| + i|l,l)(l,l|. (47) 

Already at this point, we conclude this state is obviously 
entangled. 

Let us consider the remaining operator R = g — a, 

i?=i(|0,0)(l,l| + |l,l)<0,0|). (48) 

Obviously it holds £ = 2R, see Eq. ©. Again the solu- 
tions of the separability eigenvalue problem of R are 



±i, -L(|0) + exp(#)|l}) ® -L(|0) ± exp(-^)|l}) 



These are states which fulfill Eq. flU} or Eq. l|45|) . Since 
only 2 independent bases are needed, let us choose 
exp(z0) = 1, i, —1, —i. We obtain states \s m ,S n ) together 
with the following order, 

for g = +±: |So,So>. |Sl,S 3 >, |*2,-62>, |S3,«1>) 



for g 



\So,S2), Si, Si), 02, So), |S3,S3) 
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G =4 



(49) 



(50) 



The resulting problem is given by G, g, 

9= \(h 1,1, 1,-1, -1,-1,-1) 

/ 4 1 1 1 1 \ 
14101010 
14 10 10 1 
10 14 10 10 
10 14 10 1 
10 10 14 10 
10 10 14 1 
V 1 1 1 1 4 / 

The eigenvalues A and eigenvectors p\ are 

for A = 2: p 2 = (1, 1, 1, 1, 1, 1, 1, if 
for A = 0: p = (1, -1, 1,-1, 1, -1, 1, -I) 1 
for A = 1: p M = (1,0, 1,0, -1,0,-1, 0) T 
P!, 2 = (0,1, 0,1, 0,-1,0, -if 
p li3 = (1,0, -1,0, 0,0,0, of 
pi >4 = (0,1, 0,-1, 0,0,0, of 
pi, 5 = (0,0, 0,0, 1,0,-1, of 
= (0,0, 0,0, 0,1,0, -if . 



The vector of separability eigenvalues can be expanded 



as g = jpi.i + jPi,2- Thus, the problem Gp = g can be 
solved for example for 



1 - 1 

V 
1 



P =jPl,l + jPl,2 + ^P0 



(1,0,1, 0,0,-1, 0,-l) J 



Therefore the operator R is given as 



(51) 



2R =<£ = |So,S )(s ,So| + |S2,S2)(s 2 ,S 2 | 

- |Si,Si)(Si,Si| - |s 3 ,S 3 )(s 3 ,S 3 |, (52) 

which is the form as derived in Sec. IIIII Any additional 
term proportional to pb does not decrease the number of 
vectors needed to generate A. Thus the minimal number 
of vectors needed to generate the state g - together with 
1 0,0) and |1,1) - is six. 
Now we can calculate the optimized solution as 

- pEp -. 1 -. 1 -. 
PoPa 4 4 

= J(1, 1,1, 1,-1, -1,-1,-1). (53) 

Thus an optimal decomposition of the state g = 
|$+)($+ 1 in terms of separable quantum states is 



J 



g=\<P + )(<S> 
1 



(|0,0)(0,0| + |0,0)(0,0|) 



(|So,S )(So,S | + \S 1 ,5 3 )(S 1 ,S 3 \ + |S 2 ,S 2 )(S2,S2| + |s 3 , Si) (s 3 , 5i | ) 
(|5o,S 2 )(So,S 2 | + |Si,Si)(Si,Si| + |s 2 ,S )(s 2 ,So| + |S3,S3)(S 3 ,S 3 |) 



(54) 



r 



This optimal decomposition of the Bell state in terms of 
separable ones contains negativities. These negativities 
are directly related to the entanglement of the state under 
study. It differs from the first decomposition given in 
Eq. (dSD. 



V. SUMMARY AND CONCLUSIONS 

Starting with the fundamental superposition principle 
of quantum theory, we have demonstrated its equiva- 
lence to the concept of quasi-probabilities. The repre- 
sentation of entangled states by quasi-probabilities, PEnt, 
and separable quantum states has been considered and 
compared with the Glaubcr-Sudarshan representation of 
nonclassical states in terms of classical (coherent) ones. 



We have discussed aspects of negativities in the quasi- 
distributions PEnt in relation to nonclassical Glauber- 
Sudarshan distributions. Especially, we have explained 
how entanglement is embedded into the class of all non- 
classical phenomena for the case of a two-mode harmonic 
oscillator. 

Based on the general structure of composite quantum 
states, we have studied the reconstruction of the quasi- 
probability distributions PEnt of any quantum state from 
experimental data. Starting from a Bell state, we have 
shown that for all further considerations only this exam- 
ple needs to be understood. The reconstruction proce- 
dure is based on the spectral decomposition of the statis- 
tical operator and the Schmidt decomposition of each 
eigenvector. Entangled states necessarily have quasi- 
propability distributions with negativities, independent 



9 



of the decomposition of the state. The representation of 
a quantum state in terms of separable ones and quasi- 
distributions is not unique. A positive PEnt proves sepa- 
rability, but negativities are insufficient to verify entan- 
glement. 

To obtain a unique description by quasi-probabilities, 
we have introduced an optimization procedure which de- 
livers an optimal quasi-distribution PEnt, opt- The idea is 
based on a minimal amount of negativities in the result- 
ing quasi-distribution. We have derived our optimiza- 
tion on the basis of the recently introduced separabil- 
ity eigenvalue equations. Negativities in the resulting 



quasi-probability, PEnt, opt, are sufficient to verify entan- 
glement. The optimization delivers a positive PEnt, opt 
for separable quantum states. This approach has been 
applied to two fundamental examples. 
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